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We consider a variational problem for the two-dimensional (2D) Heisenberg and XY 
models, using a trial state which is constructed as a 2D product of local weights. Variational 
energy is calculated by use of the the corner transfer matrix renormalization group (CTMRG) 
method, and its upper bound is surveyed. The variational approach is a way of applying the 
density matrix renormalization group method (DMRG) to infinite size 2D quantum systems. 

§1. Introduction 

As a precise numerical method for one-dimensional (ID) quantum systems, the 
density matrix renormalization group ©MRGT) has been widely applied to unsolved 
problems in condensed matter physics. EP'EP'Ep'q) If the linear dimension of the system 
is not so large (~ of the order of ten), the method is also applicable to 2D quantum 
models by use of a mapping from a finite size 2D cluster to a ID chain that contains 
long-range interactions. u 'EP 

On the other hand, it is difficult to apply DMRG to infinitely large 2D quantum 
systems. Apparently the mapping from 2D lattice to ID chain is inapplicable if the 
system size is infinite. In addition, the decay of the density matrix eigenvalue is 
very slow in higher dimensions, LP' tf and the phenomenon prevents to obtain a good 
renormalization group transformation that maps half (or quarter) infinite 2D system 
into a block spin. A naive extension of DMRG formulation to higher dimension 
shall encounter very poor numerical results, as was reported in an application of 
DMRG to 3D classical systems. § 

Stepping back to the DMRG formulation for ID systems, the numerical effi- 
ciency of DMRG partially comes from its variational structure, where a trial state is 
represented as a product of orthogonal matrices. Ii2rll2r 113)' EJ) It is possible to extend 
such a construction of variational state to 2D quantum and 3D classical systems. For 
these higher dimensional systems one has to prepare a trial states that are repre- 
sented as 2D product of local weights. For example, Martm-Delgad, et at employed 
the 6- vertex model as a trial state for 2D lattice spin/electron systems. EiP Okunishi 
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** Existence of a good renormalization group transformation is not excluded, in case that the 
transformation is created some method other than the diagonalization of the density matrix. 
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and Nishino considered an extension of the Kramers- Wannier approximation^ to 
the 3ELIsing model, representing the trial state as the 2D Ising model under magnetic 
field, £3) or more generally as the interaction round a face (IRF) model. E^P Let us call 
these variational approaches as 'the tensor product variational approach' (TPVA) in 
the following. 

In this paper we investigate the numerical efficiency of TPVA when it is applied 
to the square lattice 5 = 1/2 Heisenberg and XY models. We employ an isotropic 
and uniform IRF model as the trial state, that contains 3 adjustable parameters. 
Our trial state includes the one parameter variational state proposed by Suzuki and 
Miyashita for the study of square lattice XY model. @).@) The local construction 
of the variational state enables us to numerically calculate the energy expectation 
value by use of DMRGEP or by Monte Carlo simulations. c3 

In the next section we explain the construction of the trial state, and the way of 
calculating the energy expectation value. We show the numerical result in §3, and 
discuss a way of improving the variational formulation in the last section. 

§2. Construction of the Variational State 

We consider the 5 = 1/2 XXZ model on the square lattice as an example of 2D 
quantum systems. Its Hamiltonian is represented as 

H = ( s Z S r' + S^S^, + aS*S*,) , (2-1) 

{rr'} 

where {rr'} denotes pairs of neighboring lattice sites. The parameter a is chosen 
to be either 1 (the Heisenberg model) or zero (the XY model). Since the lattice is 
bipartite, the following Hamiltonian 

H = (SrS*, - S?S* + aS*S*,) , (2-2) 

{rr'} 

has the same same energy spectrum as H. In the following we treat H instead of H 
for the purpose of simplifying the variational formulation.u Our interest is in finding 
out a good variational function that minimizes the energy expectation value 

within the restriction where \W) is represented as a product of local weights. 

Let us introduce a notation a r = 2S£ = ±1, and write the trial wave function 
as ^({ct}) = where {a} represents a spin configuration of all the spins on 

the 2D lattice. We employ a trial wave function in the form of the uniform product 
of local weights 

#({a}) = Y[w( J r+ j + J j . (2-4) 



* One has to prepare two different local weights for H and multiply them alternatively when 
constructing the variational state. Note that the mapping from H to H has nothing to do with the 
elimination of negative sign from the local weight. 
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The weight W is a function of 4 neighboring spin variables a v , c r _j_j, o" r+ j, and o" r+ j + j 

where i and j represent the unit lattice vector to X and Y directions, respectively. 
According to the symmetries of the Hamiltonian H, the local weight has only 3 
independent parameters, that are 



1 = W 
a = W 
b = W 
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(2-5) 



where we have written up and down spins by '+' and '— ', respectively. Since we have 
constructed the trial state as a uniform product, the minimization of the variational 
ratio (Eq.(2-3)) is equivalent to the minimization of the local energy for an arbitrary 
bond 



A(a, b, c) 



(2-6) 



where r and r' are neighboring lattice points. In addition, the product structure of 
the trial state enables us to obtain the denominator 



<w = n 
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(2-7) 



as the partition function of the isotropic IRF model, which is specified by the local 
Boltzmann weig ht W 2 . Thus the norm of the triaLstate can be accurately calculated 
by use of DMRG applied to 2D classical systems. c3 Similarly the numerator is also 
a partition function of an IRF model that has additional structure around the bond 
{rr'} in Eq.(2-6), and can becalculated with sufficienLnumerical precision by DMRG 
as proposed by Hieida et al© We use CTMRG,©'® which is a variant of DMRG, 
for the calculations of the following result. 



§3. Calculated Result 



The minimum of the variational energy A(a, b, c) in Eq.(2-6) can be detected by 
way of the parameter sweep for a, b, and c. Table I shows the optimal parameter 
sets that give the lowest variational energy A min for the isotropic Heisenberg model 
(a = 1) and the XY model (a = 0). The diagonal spin correlation S^S^., and the off- 
diagonal one S^S*/ between the neighboring sites are also shown. In both cases the 
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Table I. Variational Energy and the value of optimal variational parameters. 





4A min 


nx nx 




a 


b 


c 


Heisenberg 


- 1.3089 


- 0.1245 


- 0.0782 


0.8342 


0.7483 


0.5075 


XY 


- 1.0848 


- 0.1356 


- 0.0383 


0.9420 


0.8796 


0.6936 



variational state is disordered, since there is no phase transition from the parameter 
limit a = b = c = lto the shown parameter cases. Even at the Heisenberg point, the 
state is not isotropic, as is observed from the difference between S^S^i and S^S*,. 



Table II. Energy expectation values. 





Heisenberg 


XY 


SM0 
TP\44 


- 1.3089 

- 1.33888 


- 1.0743 

- 1.0848 

- 1.09764 



Let us compare the calculated variational energy with two representative ground 
state energy estimations. Table II shows our result, a recent Monte Carlo (MC) 
result bySandvik, and the variational energy calculated by Suzuki and Miyashita 
in 1978.113 Our result for the Heisenberg (or the XY) model is 2.3 % (or 1.2 %) 
higher than the MC result. 

§4. Discussions 

As an application of DMRG for infinite size 2D quantum systems, we employ 
TPVA for the square lattice S = 1/2 Heisenberg and XY models. The obtained vari- 
ational state is disordered even for the Heiseberg point, the result which suggests the 
loss of the antiferromagnetic correlation. The problem may be improved by increas- 
ing the variational parameter, as Nishino et al have introduced auxiliary variables 
to the local weight when they applied TPVA to 3D classical models. El' EP'EI 
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